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Abstract 

In a toy model with gases of membranes and strings wrapping over a two-dimensional internal 
torus, we study the stabilization problem for the shape modulus. It is observed that winding modes 
of partially wrapped strings and momentum modes give rise to stress in the energy momentum 
tensor. We show that this stress dynamically stabilizes the shape modulus of the two torus. 
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I. INTRODUCTION 



In usual Kaluza-Klein picture, the extra dimensions predicted by string/M theory are as- 
sumed to be compact and small. The four dimensional low energy content of the theory (like 
the spectrum of light particles, the strengths of the coupling constants etc.) is determined by 
the internal space. Although phenomenologically viable models can be obtained in Calabi- 
Yau or G2 manifold compactifications, in its current formulation of string/M theory there 
is no principle which would single out a specific compact space. This is a major problem 
which prevents one to test the validity of string/M theory by low energy phenomena. 

Besides the vacuum selection problem, one should also confront the issue of the cosmolog- 
ical evolution of the extra dimensions. It is well known that there is a limit to the change in 
the size of the internal space following the epoch of big bang nucleosynthesis [1]. Moreover, 
the size and the shape of the compact space are parametrized by moduli scalar fields which 
are in general dynamical. Therefore, to avoid any conflict with observations, moduli fields 
should somehow be fixed in string/M theory. 

Brane gas cosmology, which was initiated in [2] and then developed in [3], may offer a 
solution to the stabilization problem in a cosmological context. 1 In [26-28], it is found that 
string winding and momentum modes can stabilize the radii of an internal torus 2 . Moreover, 
as shown in [29, 30], this mechanism survives in the presence of linearized cosmological 
perturbations. In [32, 33] we generalize these results to higher dimensional branes and 
find that the volume modulus of an internal Ricci flat manifold can be stabilized by a gas 
of branes wrapping over it. This extension is important for two reasons. Firstly winding 
strings do not exist in the spectrum when the first homology class of the internal space is 
zero and secondly they may not be able stabilize volume modulus when the topology is not 
equal to the product of 1-cycles. In [34], we also show that long- wavelength cosmological 
perturbations does not alter this result. 

1 Different aspects of brane gas cosmology are studied in the literature: T-duality invariance, a possible 
loitering phase, generalization to curved Ricci flat manifolds and M-theory are discussed in the papers 
[4], [5], [6] and [7], respectively. Cosmological solutions with brane gas sources are constructed in [8- 
14]. Annihilation of winding branes which is crucial for the Brandenberger and Vafa decompactification 
mechanism is studied in [15-20]. Inflation in the context of brane gas cosmology is considered in [21, 22] 
and the stabilization of extra dimensions in a democratic wrapping scheme is demonstrated in [23-25]. 

2 See [35] for some concerns about dilaton stabilization by string gases and [36] for an alternative approach 
to stabilization. 
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In a recent paper [37], the stabilization of extra dimensions is reconsidered in the presence 
of a string gas carrying two- form flux and it is found that, in addition to volume modulus, 
the flux and the shape moduli are also dynamically stabilized. In this paper, we address 
the stabilization problem for the shape modulus without introducing any flux. Following 
our earlier work [32] , we consider a toy model in Einstein gravity with membrane and string 
gases wrapping over an internal two-torus, but this time we take shape to be dynamical. 
In that case, we observe that winding modes of partially wrapped branes (strings in this 
model) and momentum modes give rise to stress in the energy momentum tensor. We will 
show that this stress is capable of dynamically stabilizing the shape modulus. 

As for the size moduli, we will see that when the observed space is three dimensional the 
winding and the momentum modes stabilize them as in [32]. When the dimension is greater 
than three, there appears a complication which was pointed out earlier in [9], i.e. partially 
wrapped branes force the transverse compact directions to expand. For stabilization of the 
size moduli in that case, one should restrict the relevant physical parameters so that this 
does not produce an instability. (In [38], which appeared while our work was being finalized, 
it is also found that special massless string modes can stabilize volume, shape and dilaton 
moduli.) 

The organization of the paper is as follows. In the next section we calculate the energy 
momentum tensors and obtain a background solution corresponding to a fixed rectangular 
torus and an expanding observed space. In section III, we study the linearized fluctuations 
around the background geometry and find that the shape modulus is dynamically stabilized. 
In this section, we also discuss the stabilization of the size moduli. In section IV, we conclude 
with brief remarks and future directions. 

II. ENERGY MOMENTUM TENSORS AND THE BACKGROUND SOLUTION 

We consider an (m + 3)-dimensional space-time with the metric 

ds 2 = -dt 2 + e 2B dx t dx i + G ab dy a dy\ (1) 

where i,j = l,..,m and a, b = 1,2. The m-dimensional observed space is spanned by x l 
and y a = (y 1 ^ 2 ) are coordinates on an internal two-torus T 2 so that y a ~ y a + 2n. In a 
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cosmological setup one should take 



B = B(t), G ab = G ab (t). 

A convenient parametrization of the metric on T 2 is 

Rj R X R 2 sin 9 



(2) 



G 



ab 



R X R 2 sin# R 2 2 



(3) 



(4) 



where R\ and R 2 are the radii of y 1 and y 2 , respectively, and 9 is the shape modulus. Note 
that 9 = corresponds to a rectangular torus. We choose 

R 1 cos(9/2) R 2 sm(9/2) 
R 1 sm(9/2) R 2 cos(9/2) 

as an orthonormal basis in the tangent space (the hatted and unhatted indices refer to the 
tangent and coordinate spaces, respectively). 

Up to a certain constant prefactor, the energy of a membrane wrapping over T 2 is equal 

to 

E w = Vdet G. (5) 

In the observed space the winding energy is concentrated at the position of the membrane 
and this gives a delta function singularity. But for a gas of such branes the delta function 
is smoothed out in the continuum fluid approximation [9]. The energy spectrum of the 
momentum modes corresponding to small vibrations of a membrane on T 2 is equal to the 
eigenvalues of the Laplacian of the metric G ab which are characterized by two integers 
n a = (ni,n 2 ). The energy of a mode described by the wave-function e tyaUa is given by 



E m = ^G ab n a n b . 



(6) 



Note that E m is equal for two modes having quantum numbers n a and — n a . 
For strings wrapping over yi and y 2 , the winding energy is equal to 



E s — \J Gn + \J G 22 . 



(7) 



As in the membrane case, the energy is distributed smoothly for a gas in the continuum 
approximation. String momentum modes, 3 which are labeled by n a = (nx,n 2 ) with ri\ = 



3 Here we consider semiclassical quantization of strings in a (non-covariant) physical gauge. Thus the 
tachyon is absent in the spectrum and there appears no critical dimension. 
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or ri2 = 0, can be viewed as a subset of membrane momentum modes and their energy is 
still given by (6). The corresponding wave-functions are e iylni and e iy2ri2 , which represent 
small vibrations of strings wrapping over y l and y 2 , respectively. 
From (5), (6) and (7) the energy densities can be calculated as 

E 



P 



Vole' 



(8) 



where Vols = e mB Vdet G is the total spatial volume. It is well known that (see e.g. [39]) 
one can write a matter Lagrangian 

C m = -2p 



for an energy density p. Coupling C m to Einstein-Hilbert action 

S = J y/=j[R + £ m ], 

one obtains the field equations 



where 



T 

-L II. 



d 



-g£ r . 



(9) 



(10) 



(11) 



(12) 



is the energy momentum tensor. 

Using (5) in the above formulas, the energy momentum tensor for membrane winding 
modes can be found as 



N w e 



-mB 



Too 

T = 

T N p-mBn 

± ah ly w^ ^abi 



(13) 



where N w is a numerical factor containing the tension and the number density of membranes. 
For string winding modes (7) implies 



J 00 — e 
Tij = 0, 
T ab = -e~ mB 



v/deTG s v/deTG 



, r GlaGib G 2a G2b 

IV Sl _— , ; = + iV s2 - 



(14) 
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where N s i and N s2 are proportional to the tension and the number densities of strings 
wrapping over y 1 and y 2 , respectively. Finally, the energy momentum tensor for a gas of 
string and membrane momentum modes with fixed n a can be obtained from (6) as 



Tn = o, (15) 
Th — N p~ mB naUb 

where N m is the number density Note that (15) obeys g^ u T^ u = and moreover = 0, 
therefore momentum modes behave like a gas of massless particles confined in the compact 
space. One can verify that all tensors (13), (14) and (15) are conserved 

V^T"" = 0, (16) 

so that the field equations (11) are consistent. 

It is clear from (5) and (7) that winding modes become heavy as T 2 expands and in 
thermal equilibrium they are expected to decay into lighter excitations. However, in usual 
brane gas scenario winding modes are assumed to fall out of thermal equilibrium, so one can 
take N w in (13) and N sl , N s2 in (14) as constant numbers. Similarly, from (6) the momentum 
modes become heavy as T 2 shrinks. Thus, for momentum modes to stop the contraction 
of the internal space, there should be a mechanism for them to survive the annihilation 
process. For a string gas (in dilaton gravity) the annihilation rates of the winding and 
the momentum modes are equal to each other by T-duality invariance of the field and the 
Boltzmann equations [19]. Therefore T-duality dictates that momentum modes should also 
fall out of thermal equilibrium like winding modes. Although the toy model we consider 
does not have T-duality invariance, the final scenario (obtained by including all possible 
excitations) is expected to have this symmetry [4]. Viewing the model here as a part of the 
complete picture, we assume that momentum modes also fall out of thermal equilibrium and 
thus N m in (15) is a constant number density. 

It is possible to give an alternative argument without invoking T-duality invariance that 
momentum modes should exist in the spectrum when the internal space contracts too much. 
Branes are dynamical objects and one would expect them to fluctuate when they are forced 
to collapse totally. This is crucial for stability since otherwise any unwrapped brane propa- 
gating in the bulk would undergo a complete collapse under the influence of its own tension. 
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For momentum and (partially wrapped) string winding modes, (14) and (15) imply non- 
vanishing stress 4 in the energy momentum tensor since in general ^ when a ^ b. It 
is important to refer to an orthonormal basis in determining the existence of stress. For in- 
stance, the off diagonal components of the winding energy momentum tensor of membranes 
(13) is non- vanishing in the coordinate basis. As we will see, however, (totally wrapped) 
membrane winding modes do not play a role in the dynamics of the shape modulus. 

To obtain the background solution we assume 

B = B(t), G ab = const. (17) 

Eq. (11) then yields 



(m + 1) Y v / detG 



e -mB 



(m-l)N sl ^_ mB {m-l)N s2 JG^ 2 ^_ mB 



(m+1) v/detG (m + 1) Vd^tG 

B \ m B 2 = 3Nw c~ mB I 2Nsl ^l]- e ~mB + 2N s2 \/G^ e -mB QgX 

m + 1 (m + 1) v/delG (m + 1) y^tG 

(m — 2)N W G la Gi b G 2a G 2 b 

■ ^ab + iV sl — — + iV s2 



(m+1) ^GnVdrtG ^/G^VdrtG 

n\ b 'n 



2N sl VG^ 2N s2 ^G^ 



(m+ljVdetG (m+l)VdetG 



^ + E ^ ; ifih , (20) 



where N$ is the number density for the i'th momentum mode with quantum numbers n$. 

For a rectangular torus with 9 = (i.e. G± 2 = 0) the energies of the modes with (ni, n 2 ) 
and (ni, — n 2 ) are equal to each other. Therefore, these modes should have the same number 
density. This implies that the off-diagonal component of (20) is identically satisfied after 
the sum over the pairs (ni,n 2 ) and (ni, —n 2 ). Taking moreover Ri = R 2 = R , the number 
densities of strings should also be the same, i.e. N s i = N s2 = N s . Eq. (20) then gives 

n^0^^ + 2j^^^ N{ ,^i^ (21) 

Note that the cubic equation for Rq has one unique positive root for m > 2, and for m < 2 
the ansatz is not valid. From now on, we restrict the dimension of the observed space m > 2. 



4 This stress was not observed in [32] since only the volume modulus was taken to be dynamical there. 
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Eq. (19) can be solved to get £ 



B(t) = — \n(bt), (22) 
m 



where 

b2 = 3m N w 2m N s 

2(m+l) (m+l)Ro' { ' 

One can now check that (18) is identically satisfied. Therefore the metric 

ds 2 = -dt 2 + {bt^d^dx 1 + Rl dy a dy a (24) 

solves all field equations and describes a background supported by brane winding and mo- 
mentum modes. The moduli fields of this solution are R±, R 2 and 9, corresponding to the 
radii of y 1 , y 2 and the angle between them. They have the vacuum values R\ = R 2 = Ro 
and 9 = 0. From the energy-momentum tensors (13)-(15), we see that the winding and the 
momentum modes apply negative and positive pressures along the compact directions, and 
they behave like pressureless dust in the observed space. Therefore, it is not surprising that 
in (24) the observed space expands exactly with the same power for pressureless dust and 
the internal directions are stabilized under the action of positive and negative pressures. 

III. STABILIZATION OF THE MODULI 

Having obtained the background solution, we now study linearized perturbations around 
it. The fluctuations are denoted by SB, SR±, 5R 2 and 59, which are taken to be the functions 
of time. We carry the calculation in the orthonormal basis (4). Variation of (13), which is 
the winding energy momentum tensor of membranes, gives 

5T 66 = -m N w e~ mB 5B, = mN w e~ mB 5B 5 ah . (25) 

From (14), for string winding modes we find 

5T 66 = e- mB 5B-^ 2 e- mB [5R 1 + 5R 2 ] , 

Wah = ^ e~ mB ^ 5 ab + § e~ mB [SRi 5 a2 5 b2 + 5R 2 5 al 8 bl ] 

_^± e -mB [ 6al6b2 + 6a26bl] ft. (26) 

Kq 



5 By defining a new time coordinate it is possible to obtain the most general solution of (19), which 
asymptotically becomes (22). Therefore, (22) describes the late time behavior we are interested in. 
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On the other hand, the perturbation of (15) for momentum modes yields 



\n\N„ 



oo 



-m 



-mB 



5B 



N, 



in ( mB 



n a n b Se a b 



\n\K 



Rq 



5T : 



ab 



Rt\n\ 
2N, 



n a n h e~ mB [m R 5B + 5e c c ] + 



Ro\n\ 



5e\ 

n a n b e~ mB n c n d 5e c d 



-^-^n c n (b 5e c a) , 



(27) 



(28) 



where 



5e a h 



(29) 



5Ri Rode/2 
R 56/2 5R 2 

\n\ = \Jn\ + 71% and all operations on indices are performed with the Kronecker delta func- 
tion. Note that (27) and (28) are invariant under n a — > — n a . 

The number densities of momentum modes are expected to be distributed according 
to Boltzmann weight. This does not contradict with the assumption that they fall out 
of thermal equilibrium, since just before decoupling they are expected to obey Boltzmann 
distribution and then evolve accordingly. Therefore it is sufficient to consider the modes 
having minimum energy, since the contributions of others are exponentially suppressed. In 
our case the minimum energy excitations have quantum numbers n a = (1,0), n a = (—1,0), 
n a = (0, 1) and n a = (0,-1), which should have the same density N m since they have the 
same energy. Note that both strings and membranes give rise to such momentum modes 
and N m denotes the total number density. 

Neglecting all other momentum modes, the common radius of the internal torus obeys 

<m ~ 2)N ^ + im - 3)N 'H^2N m . (30) 



(m + 1) 

Rewriting the field equations in the form 



(m+ 1) 



R 



I IV 



- )1V 



^ (m + iy 

and using (25)- (28), we find the following second order equations 



4 • 2 ■ 
SB + - SB + — SY 
t mt 



m(3N w + AN 8 ] 
(m + l)R (bty 



SY + -SY 



SZ + -SZ = 

t 



so + -se = 



[m 



S)N S 6iV„ 



1)R 







R 3 



N s | 2N m 

Ro Ro 

2N S AN m 

Rri 



1 



R* 



(btf 
1 

W 2 



(bty 
sz, 



5B- 



SY, 



2N X 



[m+l)Ro (bt)' 



SY, 



59, 



(31) 

(32) 
(33) 
(34) 
(35) 
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and the first order initial constraint 



( m _ i)SB + 5Y + SB + fcLJl 5 r = o, (36) 

t mt 



where 



5Y = i- + fiRa] , (37) 

5Z = i- - <5i? 2 ] . (38) 
-n-o 

As a consistency check of these equations we verify that the time derivative of (36) is 
identically satisfied upon using the field equations (32) and (33). Note that 5Y can be 
solved from (33) which in turn can be used in (32) to fix SB. The perturbations 5Y and 5Z 
are related to volume R1R2 and fraction R1/R2 moduli. 
Let us first consider the equation (35) for 59. Writing 

59 = t k (39) 

one finds that k obeys 

k 2 + k + c = 0, (40) 
where c = 2N s /{b 2 R ) + AN m /{b 2 Rl). The roots are 

k = - l -±i^ (41) 

where 

A = ^1. (42) 

Using the value of b in (23), we see that A is a positive reel number for m > 2. This gives 
two solutions 

cos(Alnt) sin(Alnt) , . 

» = c^-^ + c^-^, (43) 

which shows that 59 — > and the shape modulus is stabilized. 

For 5Y and 5Z, one can repeat the same analysis and get equation (40) for the corre- 
sponding power k, where the constant c is the negative of the total number multiplying 5Y/t 2 
or 5Z/t 2 in the right hand side of the linearized equation. It is necessary for stabilization 
that Re k > and by (40) this happens when 

c > 0. (44) 
10 



From (34), we see that c = 2N s /(b 2 Ro) + AN m /(b 2 Rl) > for SZ and thus the fraction 
modulus R1/R2 is stabilized. Eq. (33) shows that the volume modulus R\R 2 is stabilized 
when m = 3 since the corresponding constant c = 6N m /(b 2 Rl) > 0. To stabilize R\R 2 for 
m > 3 one should impose that N m > R^{m — 3)iV s /(6m + 6). This condition arises due to an 
effect which was pointed out in [9], i.e. (depending on the dimension of the observed space) 
partially wrapped branes force the transverse compact dimensions to expand. The condition 
on the number densities avoids the presence of an instability. Up to this complication, which 
occurs when m > 3, all moduli fields in this toy model are stabilized by brane winding and 
momentum modes. Ignoring the string gas, i.e. setting N s = 0, we find that all fields are 
still stabilized by membrane winding and momentum modes. Neglecting, on the other hand, 
the membrane gas by setting N w = 0, SY and thus the volume modulus becomes unstable. 
Focusing on the shape modulus, we see that either the momentum or the winding stress is 
capable of dynamically stabilizing SO. 

Therefore we see that the negative terms which appear in the right hand sides of (33)- 
(35) are responsible for the stabilization of the moduli fields. It seems that the membrane 
winding modes characterized by N w does not play a role in the stabilization since N w does 
not appear in the right hand sides. This is only true at the linearized level for SRi and SR 2 , 
and the membrane winding modes help for the stabilization of the radii nonlinearly [32]. 
However the right hand side of (35), which is (/t, v) = (1,2) component of (31), is equal to 
the off-diagonal component of the energy momentum tensor in the orthonormal basis Tjg. 
From (13), we see that Tj§ = for membrane winding modes exactly. Therefore, winding 
modes of totally wrapped branes do not play a role in the dynamics of the shape modulus. 
The situation is different for partially wrapped branes. From (33)-(35) we see that winding 
string modes, which are characterized by N s , work for the stabilization of SO and SZ, but 
they try to destabilize SY when m > 3. In the linearized approximation, the contributions 
of momentum modes n a = (1, 0), n a — (—1, 0), n a — (0, 1) and n a = (0, —1) are negative in 
the right hand sides of (33)- (35), thus they try to stabilize all moduli fields. 

Although the higher level momentum modes are exponentially suppressed, one can de- 
termine their effect on the stabilization of the shape modulus. From (28), we find that 

AT AT K + ^) SB 

5Ti2= - Nm ^mrw (45) 

for the sum of the modes (ni,n 2 ) and (ni, —n 2 ). This term should be added to the right 
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hand side of (35) and since it is negative we see that the pairs of the modes (711,712) and 
(tii, — n 2 ) force 59 to be fixed. 

IV. CONCLUSIONS 

In this work, we study the stabilization problem for the shape modulus in a toy model in 
Einstein gravity with gases of membranes and strings wrapping over an internal two-torus. 
Based on the earlier results in the literature (for strings [26-28] and for higher dimensional 
branes [32, 33]), it is not surprising to see that the size moduli are stabilized. Our new 
observation is that the winding modes of partially wrapped branes and the momentum 
modes give rise to stress which is capable of dynamically stabilizing the shape modulus. We 
also notice that winding modes of totally wrapped branes do not play a role in the dynamics 
of the shape modulus. Of course the model we consider is incomplete in many ways but it 
nicely illustrates a generic dynamical behavior. Although for a two-torus there is only one 
shape field to worry about, the mechanism can be generalized to work for higher dimensional 
tori with more than one shape moduli . 

In this paper we carry the calculations in the framework of Einstein gravity. It would be 
interesting to generalize these findings to 10-dimensional dilaton gravity. In that case, the 
results would depend on the brane type. In [33], consistent field equations are obtained for 
a D-brane gas coupled to dilaton gravity and a natural extension of this work is to consider 
a dynamical shape in that toy model. 

It is known that in string theory compactifications with fluxes, all of the complex structure 
moduli including the shape of the internal manifold can be stabilized. However, it turns out 
to be difficult to stabilize the volume modulus in this framework. On the other hand, volume 
stabilization can be achieved easily in brane gas cosmology. This motivates the inclusion 
of brane gases carrying fluxes and as shown in [37] all moduli fields can be stabilized in a 
simple toroidal compactification with a string gas and ambient flux. In this work we show 
that the shape moduli can be stabilized without introducing fluxes. However, one may need 
to introduce fluxes for other moduli (like the ones related to form- fields generating fluxes), 
therefore it would be interesting to develop further the framework of brane gases with fluxes. 

We study the stabilization problem in a linearized setting. It would be important to 
extend the results to the non-linear regime without any apprxoximations. One difficulty in 
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this analysis is to determine the relative number densities of the momentum modes. For 
instance, in a square two-torus with 9 = 0, n a = (1, 0) mode is lighter than n a = (1, 1) 
mode. However, as 9 increases, the energy of n a = (1, 1) mode decreases. When sin# = 1/2, 
the energies of the modes n a = (1,0) and n a = (1, 1) become equal. As 9 increases more, 
n a = (1, 1) mode becomes the lighter one. (In [40] and [41], this relation between the shape 
of the torus and the energy of the momentum modes was used to weaken the experimental 
bounds on theories with large extra dimensions). In a non-linear analysis, fluctuations on 
9 can be large and it is difficult to determine the light species which would dominate the 
dynamics. 

The torus is a special manifold where the eigenvalues of the Laplacian are known explicitly. 
It is worth trying to extend the arguments where this information is no longer available. 
In the realistic case, to any topologically nontrivial cycle, one anticipates the existence of a 
brane gas wrapping over that cycle. The results in brane gas cosmology suggest that in such 
a scenario all moduli fields should be stabilized by brane winding and momentum modes. 
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